Controlled generation of coherent matter-currents using a periodic driving field 
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We study the effect of a strong, oscillating driving field on the dynamics of ultracold bosons held 
in an optical lattice. Modeling the system as a Bose-Hubbard model, we show how the driving 
field can be used to produce and maintain a coherent atomic current by controlling the phase of 
the intersite tunneling processes. We investigate both the stroboscopic and time-averaged behavior 
using Floquet theory, and demonstrate that this procedure provides a stable and precise method of 
controlling coherent quantum systems. 
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Introduction - Recent experimental advances in 
the creation of Bose-Einstein condensates (BECs) from 
ultracold atomic gases have stimulated huge interest 
in investigating the coherent many-body dynamics of 
trapped bosons. By superposing counter-propagating 
laser beams, it is possible to impose extremely well- 
controlled lattice potentials on these systems. The pre- 
cision and flexibility afforded by this both suggests their 
use as "quantum matter simulators" [1,] for systems of 
interest from other areas of physics such as the integer 
quantum Hall effect 0, and also permits the clean ob- 
servation of coherent lattice phenomena such as Bloch 
oscillations Q , the formation of repulsively-bound pairs 
0| , and the Mott transition 

As well as their purely theoretical interest, these sys- 
tems are also highly attractive candidates for applications 
such as quantum information processing due to their long 
coherence times. A powerful tool to control their dynam- 
ics is provided by the effect termed "coherent destruction 
of tunneling" Q , in which driving the system with an os- 
cillating field has the effect of renormalizing the intersite 
tunneling amplitude. For certain parameters of the driv- 
ing field the tunneling can even be reduced to zero, and 
thus it has been proposed to use this effect to control the 
quantum phase transition between the superfluid and the 
Mott state 0,(1. 

In this Letter we show how an oscillating driving field 
can not only be used to control the amplitude of the tun- 
neling, but also its phase. This permits the generation 
of a matter-current in analogy to the current induced 
in a conducting ring threaded by a magnetic flux. The 
creation of tunneling-phases has been studied before in 
BECs using either rotating lattices or by using the 
atoms' internal degrees of freedom to mimic a fictitious 
magnetic field [lo| . The scheme we propose is extremely 
simple in comparison, requiring only the controlled shak- 
ing of the optical lattice, which has already been demon- 
strated in experiment. It also reverses the normal role 
of CDT in a novel way, since the quantum interference 
effects which produce CDT are employed here to induce 
motion, rather than to suppress it. We study the effect 
over a range of interaction strengths, and find it to be 



present from weak interactions right up to the onset of 
the Mott state. A surprising feature is that while inter- 
actions do reduce the magnitude of the current, they do 
not introduce dephasing or dissipation as seen, for ex- 
ample, in Bloch oscillations [HI], but instead render the 
current-generation more robust. 

Model - We consider a one-dimensional (ID) optical 
lattice, in which the atoms are confined to the lowest 
Bloch band. In this case the system can be described 
very accurately [12| by the Bose-Hubbard (BH) model 

#BH = [~ J a ™ a n + Hc ] + jJ2n m {n m - 1) ■ 

(m,n) m 

(1) 

Here a m (a' m ) are the standard boson destruction (cre- 
ation) operators, n m — al n a m is the number opera- 
tor, and U is the Hubbard-interaction between a pair 
of bosons occupying the same site. The tunneling am- 
plitudes J connect nearest-neighbor sites (m,n), and we 
take h — 1. We now impose a time-dependent potential 
which rises linearly across the lattice 

H(t) = H BU + K(t) sm(ujt + 0)^2 n i > ( 2 ) 

3 

where lo is the frequency of the driving field, and K(t) 
parameterizes its amplitude. Importantly we include the 
phase of the driving field, 9, as an additional control pa- 
rameter. This form of potential can be produced by peri- 
odically phase-modulating one of the laser fields provid- 
ing the optical lattice, and has already been used in cold 
atom experiments [13| to induce CDT. 

We begin by considering the case of a driving field 
of constant amplitude, K{t) = K. The Hamiltonian 
of the system, Eq. [2l is then periodic, with period 
T = 2ir/ui. Accordingly we may use the Floquet theo- 
rem to write solutions of the time-dependent Schrddinger 
equation in the form u(t) = exp(— iet)u(t), where u(t) is 
a T-periodic function termed a Floquet state, and e is 
called a quasienergy. In the high-frequency limit, where 
co ^> (J, U) is the dominant energy scale of the problem, 
perturbative approximations to the Floquet states can 
be obtained by solving the Floquet equation for just the 
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driving potential, and then including i?BH as a pertur- 
bation. In this case, following the procedure described in 
Refs. [HE!, the Floquet states are given to first-order 
by the eigenstates of an operator Ti.it) which is identical 
to i/BH) but with periodically- varying tunneling ampli- 
tudes given by 



Jit) = J 



T 



(3) 



where F[t, t) — sin(w7j' + ff)dt' is the phase accumu- 
lated over the interval (i, r), and the +/— applies to for- 
ward/backward hopping. If we now consider the system 
stroboscopically, that is, at discrete times t = nT where 
n is integer, the time-dependence of the tunneling ampli- 
tudes disappears, and the system is effectively governed 
by a static Hamiltonian 7i(0), with the Floquet states 
Uj(0) playing the role of energy eigenvectors. Simplify- 
ing Eq. [3] then reveals that the action of the driving field 
is to rcnormalize the tunneling amplitudes as 



J cff = Je ±l{K ^ cos0 JoiK/u), 



(4) 



where Jo is the zeroth Bessel function of the first kind. 
For 9 = 7r/2 (cosinusoidal driving) this result reproduces 
the familiar Bessel function renormalization [7|, l8| of the 
tunneling. An unanticipated result, however, is that the 
hopping in general acquires a non-zero phase, which is 
maximum for 9 — (sinusoidal driving). 

It may appear suprising that 9 can induce a non-trivial 
phase, as changing 9 is equivalent merely to shifting the 
time origin. In practice, however, the driving field must 
be turned on at a certain time, which thereby does pick 
out a specific value for the phase of the driving field. 
Since the system is completely coherent, the effect of 
this initial condition is not lost during the subsequent 
time-evolution, and thus the driving phase 9 can pro- 
duce different physical results. The central result of our 
work is that if K is increased from zero sufficiently slowly, 
the Floquet states of the system are able to adiabatically 
follow [16], and thereby acquire the if-dependent phase 
<j) = Kit) cos 9/uj. As a result, a given initial state can be 
transformed into a current-carrying state by slowly ramp- 
ing the driving potential from zero to the value that gives 
the desired hopping-phase. 

Results - To verify these results, we study the be- 
havior of an iV-site BH system by numerically propa- 
gating the many-particle wavefunction under the time- 
dependent Hamiltonian @ . We focus on the case of com- 
mensurate filling, where the number of bosons is equal to 
N, so that in the limit of large U a well-defined Mott state 
exists. The rapid increase in the dimension of the Hilbcrt 
space means that we could only consider systems of up 
to N = 10, but examining the results as N is increased 
reveals that the behaviour we find is quite insensitive to 
lattice size. To probe the behavior of the system, we 
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FIG. 1: Normalized momentum density p(p)/iV for the first 
Brillouin zone of an 8-site lattice holding eight bosons. For a 
perfect superfluid the momentum density is sharply peaked, 
showing the presence of long-range coherence. As U is in- 
creased the peak flattens and broadens, until for U = 16J 
the distribution is almost flat, indicating the proximity to the 
formation of a Mott insulator. 
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This quantity is remarkably size-independent [17| , allow- 
ing results from small lattice systems to be reliably ex- 
trapolated to the thermodynamic limit. It can be ob- 
served directly in experiment by time-of-flight absorption 
imaging, and conveniently indicates whether the system 
is in the superfluid or Mott-insulator regime When 
interactions are weak (U <C J), bosons are delocalized 
over the lattice in a superfluid state, and the system pos- 
sesses long-range phase coherence. Consequently pip) is 
sharply peaked, as shown in Fig. [TJ As the interac- 
tion is increased the peaks broaden and reduce in height, 
indicating that the bosons become progressively more lo- 
calized on the lattice sites. For sufficiently large values 
of U/J the atoms localize completely to form the Mott 
state, for which the momentum distribution is completely 
flat. In ID this phase transition is quite soft, giving a 
large range of U over which pip) is peaked. 



We begin by considering the case of an intermediate 
interaction strength, U = 8J. To place the system in 
the high-frequency regime we set u> — 30J, and use a 
sinusoidal driving field (0 = 0). The precise form in which 
Kit) is increased from zero is not important as long as 
it satisfies the adiabaticity constraint, and for simplicity 
we consider a linear ramp Kit) = K$t. As noted earlier, 
we will evaluate all physical quantities stroboscopically 
at times t = nT. The system is initialized in its ground- 
state, and in Fig. [5^ we plot the expectation value of 
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FIG. 2: Current induced by a linearly-ramped driving field 
K(t) = Kot, where Ko = 0.05T -1 , in an 8-site system with in- 
teraction strength U = 8 J. The driving is sinusoidal (9 = 0). 
(a) The induced current (solid black line), shows a decay- 
ing oscillatory behavior, described well by Eq. [6] (dashed-red 
line). As well as at multiples of it, marked by the vertical 
lines, zeros of the current also occur when Jo(K/uj) = due 
to CDT. (b) As above, the dashed line indicates the current 
produced by a continuously ramped field. Holding K (t) con- 
stant after a certain time (shown schematically in the inset), 
keeps the current at a constant level; the solid curve shows 
the effect of ramping the field until t = 21T, which gives the 
maximum current 7q. 



the lattice current, / = 2JIm(a m aJ i ), as a function of the 
driving amplitude. For K(t) = the current is zero due 
to the symmetry of the momentum distribution. As K is 
increased, however, the peak in p(jp) is displaced from the 
center of the Brillouin zone due to the induced hopping- 
phase <p, p(jp) — > p(p + 4>). As a result I becomes non-zero 
due to the imbalance between the left and right-moving 
momentum components. 
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FIG. 3: Momentum distribution for an 8-site system (U = 8 J) 
under sinusoidal driving with a linearly-ramped amplitude. 
To guide the eye three Brillouin zones are plotted, and dashed 
lines indicate the evolution of the central peak. The verti- 
cal lines mark the zeros of Jo. As K increases the peaks 
steadily shift in momentum due to the induced hopping phase, 
and their amplitude reduces according to the Bessel function 
Jo(K/u). At K/uj = 2.4048 the Bessel function approaches 
zero and the system becomes a Mott insulator with a flat dis- 
tribution. Increasing K further causes J e g to change sign, 
and the peaks reappear with a shift of n (see text). This 
pattern then repeats. 



Ref. [l3|. The intricate behavior of I in Fig. [2^ thus 
arises from a combination of the shifting location of the 
peaks, together with many-particle effects arising from 
the competition between J e g and U. The roughness visi- 
ble in the current arises from departures from adiabatic- 
ity in the driving. As the ramping rate is decreased, and 
so approximates adiabatic evolution more closely, this 
roughness is progressively eliminated. 



Following its initial increase, I displays a damped oscil- 
latory dependence on K/uj. To interpret this behavior, 
we show the corresponding response of the momentum 
density in Fig. [3l As expected, the initial effect of the 
ramping potential is simply to shift the locations of the 
peaks in the momentum distribution by inducing the hop- 
ping phase. As the peaks shift, however, their amplitude 
is reduced by the Bessel function renormalization of the 
hopping amplitude As K/u> — > 2.4048, the first zero 
of Jo , the effective hopping vanishes and the system thus 
makes a transition to the Mott state [3] and the momen- 
tum density flattens. As K is increased further, peaks 
reappear in p(p), but their location is discretely shifted. 
This occurs because J e g has become negative; writing it 
as J e ff = | Jeff | exp(i7r) clearly indicates that the peaks in 
the momentum density will be displaced by ir. Predicted 
in Ref.[l8j], this shift has recently been experimentally 
observed for a weakly-interacting (U ~ 0.1J) system in 



For strong interactions, the ground-state of the sys- 
tem consists approximately of the Mott state, 1 11111 . . .}, 
with a small admixture of excited states. The dominant 
dynamical processes will be nearest-neighbor tunneling 
between the Mott state and "particle-hole" states, sep- 
arated from the Mott state by an energy gap of ~ U, 
where one site is doubly-occupied and one site is empty 
(e.g. 1 12011 ...))• Including only these processes, we can 
obtain an approximate form for the induced current 

I~2JH 2 sia(Kcos6/u)J (K/u), (6) 

where \a\ 2 is the weight of the particle-hole states in 
the interacting ground-state. For large U, \a\ 2 decays as 
~ U^ 1 as the ground state converges toward the Mott- 
state. Using a as a fitting parameter, we show in Fig. 
[2k that this expression indeed provides an excellent de- 
scription of the current. Eq. [6] reveals the two distinct 
sources for the zeros of current; (i) when K/lo = nir the 
momentum density is symmetrically-peaked at the center 
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FIG. 4: The maximum current, Jo induced in an 8-site system 
depends on the phase of the driving field, 0. In all cases Jo 
is maximized for sinusoidal driving and is zero for the cosinu- 
soidal case. Solid lines plot the dependence 7sin(i('cosf)/cj), 
where 7 is a fitting parameter, and show excellent agreement 
with Eq. [5] As U is increased, Io reduces, but even for 
U — 16J the induced current is significant. 

of the Brillouin zone, and the positive and negative cur- 
rents cancel, (ii) when the Bessel function becomes zero, 
J e ff is suppressed and so the tunneling itself is quenched. 

An important consequence of the interaction is that 
larger values of U confer increased stability during the 
ramping process. For weak interactions only extremely 
slow ramping can be used, or the system will be excited 
from its instantaneous ground state and control of the co- 
herent current will be lost. When U is large, however, the 
resulting energy gap isolates the ground-state from the 
rest of the spectrum and makes the adiabatic condition 
easier to attain, thereby allowing more rapid ramping to 
be used. 

Differentiating Eq. [6] reveals that the maximum cur- 
rent, I , occurs for K/u — 1.0311. In Fig. Eb we show 
the effect of ramping the value of K(t) up to this value, 
and then keeping it fixed (see inset), which maintains 
the induced current at its final value. The magnitude of 
the current depends only on the final value of K/uj, and 
so by regulating this value any desired value of current 
within the range ±Iq can be generated. In Fig. [4] we 
show the dependence of Iq on the phase of the driving 
9 for several different values of U. All the curves show 
the s'm(K cos 6/u>) dependence expected from Eq. [6l and 
the magnitude of the current remains significant even for 
large interaction strengths near the onset of the Mott 
transition. While smaller values of U allow larger cur- 
rents to be induced, slower ramping rates must then be 
used, giving a trade-off between the two effects in exper- 
imental implementations. 



Conclusions - We have described a means of induc- 
ing a coherent atomic current by adiabatically-controlling 
the renormalization of the intersite tunneling. The pres- 
ence of interactions both stabilizes this mechanism, and 
also introduces novel strong-correlation effects. To reveal 
this effect we have employed a stroboscopic measurement 
scheme; this implies that in experiment measurements 
must be made at well-controlled intervals, and be suffi- 



ciently rapid to reflect the system's instantaneous state. 
For typical cold atom systems this would require tempo- 
ral control of the order of milliseconds, which should be 
easily achievable (l3j ]. If the measurements have insuf- 
ficient time-resolution, it would then be appropriate to 
consider the time-average of Eq. [31 which yield the re- 
sult that ( J e fj) = JJq(K/u), and thus the hopping-phase 
vanishes and the effective hopping is proportional to the 
square of the Bessel function. For a tilted lattice, it can 
be shown that the tunneling is renormalized as J^{K /u>), 
where muj is the energy difference between adjacent sites. 
Interestingly such a dependence has been recently ob- 
served in [19( for m = 1,2. While we have focused on the 
case of bosonic atoms, this mechanism should be equally 
applicable to cold fermionic atoms or electronic systems, 
provided that they possess the required coherence prop- 
erties. 
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